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ABSTRACT 

The dynamical evolution of the large (> 100 kpc), barred spiral galaxy NGC 6872 and its small companion IC 4970 in the southern 
group Pavo is investigated. We present A'-body simulations with stars and gas and 21 cm Hi observations carried out with the Australia 
Telescope Compact Array of the large-scale distribution and kinematics of atomic gas. Hi is detected toward the companion, corre- 
sponding to a gas mass of ~ 1.3 X 10' Mq. NGC 6872 contains ~ 1.4 x IO'^Mq of Hi gas, distributed in an extended rotating disk. 
Massive concentrations of gas (~ 10' Mq) are found at the tip of both tidal tails and towards the break seen in the optical northern arm 
near the companion. We detect no Hi counterpart to the X-ray trail between NGC 6872 and NGC 6876, the dominant elliptical galaxy 
in the Pavo group located ~ 8' to the southeast. At the sensitivity and the resolution of the observations, there is no sign in the overall 
Hi distribution that NGC 6876 has affected the evolution of NGC 6872. There is no evidence of ram pressure stripping either. The X- 
ray trail could be due to gravitational focusing of the hot gas in the Pavo group behind NGC 6872 as the galaxy moves supersonically 
through the hot medium. The simulations of a gravitational interaction with a small nearby companion on a low-inclination prograde 
passage are able to reproduce most of the observed features of NGC 6872, including the general morphology of the galaxy, the inner 
bar, the extent of the tidal tails and the thinness of the southern tail. 

Key words, galaxies: interaction - galaxies: ISM - ISM: kinematics and dynamics ~ galaxies: individual: NGC 6872, IC 4970, 
NGC 6876 



1. Introduction 

. Galaxies with extended tidal tails form an outstanding labora- 
' tory to study the effect of a collision on the different compo- 
nents of a galaxy because of the sensitivity of the morphology 
and kinematics of the tails to the initial distribution of the matter 
(both lum inous and dark) and to the geometry of the i nterac- 
tion (e.g.. iDubinski et al.lll999l ISpringel & White! 1 19991) . Early 
simulations clearly showed that prograde, co-planar encounters 
\ are more efficient than retrograde encounters at triggering loiig 
thin tails (e.g., Holmberg 1941, Pfleiderer & Siedentopf 19611 
iPfleidereil IT9631 Pfoomre & Toomreill972l) . The tails are made 
of stars and gas that have been thrown out from the galactic 
disks during the gravitational interaction which may lead to the 
merger of the two disks. Spectacular examples of nearby galax- 
ies at differe nt stages of an interaction can be found in the Arp 
atlas (Il966h and the Arp & Madore catalog tll987h . The tails 
are often gas-rich and their kinematics can be traced out to 
large radii by 21 cm hne observations of atomic hydrogen (e.g., 
iHibbard & van Gorkom 1996). Some tidal tails contain massive 
self-gravitating concentrations of matter which may be the pro- 
genitors of the so-called tidal dwarf galaxies (lEknegreen et al.l 
I1993L iBournaud et al.ll2003l iBournaud et al.ll2004 . Tidal tails 
develop not only in merging systems, but also in flyby encoun- 
ters, where the interacting galaxies do not form a bound system. 
It is of particular interest to study tidal tails at an early time of the 
interaction when the two galaxies are well separated and have a 
clear morphology. A few 10^ yr after closest approach, the tails 
are already well developed and it is possible to reconstruct the 
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dynamical history of the interaction and examine its influence 
on the different components of a galaxy and the level of induced 
star formation. 

The southern galaxy NGC 6872 is one of the largest spi- 
ral galaxies known. Star formation is traced all along the arms, 
which extend over more than 100 kpc at our adopted distance of 
61 Mpc. The galaxy belongs to a small group, Pavo (see Fig. la). 
NGC 6872 is likely to be affected by tidal perturbations from 
the nearby com panion IC 4970, a small lenticular galaxy located 
in to the north. iMachacek et alj ( l2005l) have discovered a more 
than 100 kpc long X-ray trail extending between NGC 6872 and 
the neighbor galaxy NGC 6876, the dominant elliptical in the 
group which lies about 8' (~ 142 kpc) to the southeast. The ra- 
dial velocity of NGC 6876 is about 800 km s ' lower than that of 
NGC 6872 (see Table 1). The X-ray trail is hotter (~ 1 keV) than 
the undisturbed Pavo intergalactic medium (~ 0.5 keV) and has 
a low metal abundance. The authors interpret the trail as partly 
due to gravitational focusing of the intracluster gas into a Bondi- 
Hoyle wake, as the spiral galaxy moves supersonically through 
the intracluster medium; they point out that the trail could also 
consist of a mixture of intracluster gas and gas removed from 
NGC 6872 by turbulent viscous stripping. 

The spectacular Very Large Telescope (VLT) multicolor im- 
age of NGC 6872/IC 4970 displays striking contrasts between 
the inner region, especially the straight northern arm, and the 
blue diffuse tidal tail to the north-east (ESO Press Release 
20b/99; see also Fig. lb, which shows the VLT blue-band im- 
age on a grey scale). A bar is clearly seen in the 2MASS J, H, K 
images d Jarre tt et al.ll2003l) . Fabry-Perot Ha observations have 
revealed the presence of ionized gas at the tip of both tails 
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Note NGC 6872 


IC 4970 


NGC 6876 


Other names 




ESO 073-IG 32 


ESO 073-IG 33 ESO 073-IG035 






VV 297a 


VV 297b 








AM 201 1-705 






Q'(J2000) 


(1) 


20M6™56'!91 


2046'"57>;87 


2ohi8"U9U5 


(5(J2000) 


(1) 


-70°46'04'.'5 


-70°44'57'.'3 


-70°51'31'.'7 


Type 


(2) 


.SBS3P. 


.LA.-P* 


SBO 


Systemic velocity (2) 


4818 km S-' 


4727 km s"' 


4010 km s"' 


Distance 


(3) 


61 Mpc 






Scale 




17.75 kpc arcmin" 
296 pc arcsec"' 


1 




D25 


(2) 


6:025=106.9 kpc 


0;676=12kpc 


2f8 


Position angle 


(2) 


66° 


6° 


79° 


Axis ratio 


(2) 


0.288=cos(73°) 


0.323=cos(7r) 0.786=cos(38°) 




(4) 


1.67 Jy 






/lOOfim 


(4) 


6.61 Jy 






^FIR 


(4) 


28 K 






^FIR 


(4) 


1.6 X lO'^Lo 






Lb 


(2) 


1.67 X 10" Lo 


1.24 X lO'^Lo 


1.07 X 10"Lo 


Hi mass 


(5) 


1.41 X lO'^Mo 


1.3 X 10' Me 




H2 mass 


(6) 


9.6 X 10** Mo 







Notes: (1) Posi tion of NGC 68 7 2: iRussell et al] jl990h . uncertainty of 2'.'5. Position of IC 4970: iLaubertsI jl982l). uncertainty of 7' 5. Position 
of NGC 6876: 'Sfautski e et Zl ([20061), uncertainty of 1'.'25. (2) Systemic velocity of NGC 6876: iDavoust & Consider^ ( Il995h : other data: 
Ide Vaucouleurs et al. ( 1991]). The blue-band luminosities have been derived from the blue-band magnitudes corrected for extinction, S^, using 
the following relation: logCLg/Lo) = 12.192 - 0.4B° + log(£)/lMpc)-, consistent with the blue magnitude for the Sun M" = 5.48. 

(3) Consistent with Hp = 75 km s"' Mpc"' . 

(4) lMoshir & et al] ( Il990h . The far-infrared luminosity Lfm, in solar luminosities, has been derived from the flux densities at 60 and 100 yum listed 
in the IRAS catalogue using the relation: Lpm = 3.94 x 10^ Z)^(2. 58/^0 + /loo) where D is the distance in Mpc, /go and /joo are the flux densities 
in Jy. The dust temperature Tpm has been derived from /go and /oo by fitting a black body and assuming that the dust emissivity is inversely 
proportional to the wavelength. 

(5) This work. 

(6) lHorelIou & BootH ( fT997l) . The H2 mass refers to the central 45" of NGC 6872 and was derived from CO(l-O) observations using a standard 
CO-H2 conversion factor. 

Table 1. Basic parameters of NGC 6872/IC 4970. 



dMihos et alj 11993]) . The location of the ionized regions coin- 
cides with that of young, blue stellar clusters seen in the optical 
images. No Ha emission was detected from the central region 
and the bar, nor from the early-type companion IC 4970. 

iBastian et alj (l2005h have studied the rich population of star 
clusters in NGC 6872 using archival VLT images in the B, V, 
and / bands, complemented with new VLT observations in the 
t/-band. They estimated the mass, age and extinction of the star 
clusters and found that most of the young massive clusters (of 
mass between 10"^ and 10^ and less than 100 Myr old) are 
located in the tidal tails or in the outer part of the galactic disk. 
The mass distribution of the star clusters follows a power-law 
with an index similar to that found in other galaxies. The authors 
also estimated the star formation rate in different regions of the 
galaxy and showed that the northern tail is forming stars at about 
twice the rate of the southern tail, and about five times the rate 
of the main body. 

iMihos et"ail (1 1993b have modeled the dynamic evolution of 
NGC 6872, as it is gravitationally perturbed during the passage 
of the small companion IC 4970 on a low -inclination, prograde 
orbit. They were able to reproduce most of the observed fea- 
tures, although the simulated tidal tails were not as thin as the 
observed ones. The model predicted an accumulation of gas to- 
ward the_£ente]vwdiich^ detected in CO emis- 
sion (iHorellou & Boothlll997h . 

In this paper, we present new maps of the large-scale 
distribution and kinematics of atomic hydrogen (Hi) in the 
NGC 6872/IC 4970 system and A^-body sim ulations with stars 
and gas. With a total Hi mass > 1.5 x 10'"Mo (iHorellou & BoothI 



IT997I mH), the pair is among the most gas-rich galaxies in 
the optically selected sampl e of ~ 60 interacting and merging 
galaxies compiled by Johan sson & Berg vall ( 1990). The single- 
dish Hi spectrum is extremely broad (~ 950 km s"' at the base 
of the Hi line) and observations at a resolution higher than the 
Parkes 15' beam were required to determine the location of the 
Hi gas. The Australia Telescope Compact Array (ATCA) in sev- 
eral configurations makes it possible to obtain data with angu- 
lar resolution (~ 40", corresponding to about ~ 12 kpc at our 
adopted distance) and a sensitivity adequate to detect extended 
Hi emission. We pointed the aiTay half-way between the large 
spiral galaxy and NGC 6876, the dominant elliptical in the Pavo 
group located ca 8' to the south-east, in order to be more sen- 
sitive to diffuse Hi gas in the Pavo group and detect traces of a 
possible gravitational interaction in the extended Hi distribution. 



We also discuss A^-body simulations with stars and gas per- 
formed using a particle-mesh code, where the gas is modeled 
as a collection of clouds that dissipate energy through inelas- 
tic collisions. Prehminary results have been presented elsew here 
(IHorellou & Koribalskill200l IHorellou & Koribalskill200l . We 
have used simulations to investigate the evolution of gas and 
stars in such a close encounter, examine the influence of the dark 
matter halo and of the geometry of the collision on the charac- 
teristics of the tidal tails and place constraints on the dynamical 
history of the observed system. 
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Fig. 1. (a) Optical image from the Digitized Sky Survey of the field around NGC 6872. North is up and east is left. iMachacek et al] 
( I2005h detected diffuse X -ray emission in the region between NGC 6872 and the massive elliptical galaxy NGC 6876. (b) Integrated 
Hi map overlaid on the blue-band VLT image. The contour levels increase from 14.4 x 10''^ cm"^ to 115.4 x 10'^ cm"^ by steps 
of 14.4 X 10''^ cm"^. The synthesized beam is shown in the bottom left corner. The frame is 155 kpc on each side for the adopted 
distance. 



2. Hi Observations 

2.1. Observations and data reduction 

Hi observations were carried out using the ATCA in five differ- 
ent configurations. Details are given in Table 2. The correlator 
provided a total bandwidth of 16 MHz divided in 512 chan- 
nels, which corresponds to a total velocity coverage of 3431 
km s"' and a velocity resolution of 6.7 km s The data were 
reduced in miriad and aips using standard procedures. After sub- 
stantial flagging (narrow-band interferences around 1400, 1403 
and 1406 MHz plus broad-band solar interference, all affect- 
ing the shortest baselines) the data were calibrated using PKS 
1934-638 which lies only 8° away from NGC 6872. The 20-cm 
continuum emission was subtracted by fitting a first order base 
level to all line-free channels. The four good Hi line data sets 
were then Fourier-transformed together using 'robust' weight- 
ing (R0BUST=0.5). To enhance the signal-to-noise ratio of the 
Hi emission we smoothed the data to a velocity resolution of 
~ 20 km s"' . All velocities are in the barycentric reference 
frame using the optical definition. The final maps obtained by 
combining the various data sets have a synthesised beam of 
42'.'4 X 38'.'8 and a noise level per channel of 0.90 mJy, which 
is very close to the theoretical value of 0.85 mJy. Moment maps 
(zeroth order, first and second order) were obtained using signals 
above 2 mJy beam"'. 

Hi masses are related to the Hi integrated intensities J Snidv 
(in Jy km s"') by: 



(1) 



M(HI)(Mo) = 2.36 X lO^D^p, J Smdv 

where Dmpc is the distance to the galaxy in Mpc. 

Hi column densities can be derived from the observed inte- 
grated intensities by using the relation: 



HI 



1 10.4 X 10^ 

ab 



Arrays, 


750D, 25/09/2001, 12h 


observing dates 


1.5D (corrupt), 17/11/2001, <6h 


and durations 


EW352, 26/12/2001, 12h 




1.5B, 07/07/2002, 12h 




750B, 03/08/2002, 12h 


Pointing position (J2000) 


2047'"38''; -70°48'47" 


Primary beam 


34:1 


Central frequency 


1399 MHz 


Number of channels 


512 


Calibrator 


PKS 1934-638 (14.91 Jy) 


Combined data: 




Synthesized beam 




- Major X minor axis 


42'.'4 X 38'.' 8 (FWHM) 


- Position angle 


+ 12° 


Channel width 


Av = 20.55 km s"' 


Noise level (Icr) 


0.9 mJy beam" 'channel ' 




1.24 X 10" cm"2 channel"' 




1.6 X 10^ Mq beam" 'channel"' 



Smdv 



(2) 



Table 2. Summary of Hi observations. 



where A^hi is the column density of hydrogen atoms in 10' cm" 
and a and b are the diameters of the synthesized beam at full 
width half maximum (FWHM) in arcseconds. 

2.2. Results 
2.2.1. Hi distribution 

Figures lb and |2] display the Hi distribution as isocontours 
superimposed on the blue-band VLT image on a grey 
scale. The atomic gas is confined to the interacting galaxies 
NGC 6872/IC 4970 and there is no evidence for the presence 
of extended Hi gas around the galaxies or in the region between 
NGC 6872 and the massive elliptical NGC 6876. We measure 
an Hi flux of 17.5 Jy km s"' toward NGC 6872/IC 4970, which 
is less than the published single-dish measurements within 



4 



Horellou & Koribalski: Stars and gas in NGC 6872/IC 4970 




Fig. 2. Details of the Hi distribution as contours overlaid on the VLT B-band image, as in the previous figure. The size of the beam 
is indicated in the lower left corner of the first two frames. Aside from the first frame, all the other frames have the same size. North 
is up and east is left. The contour levels are the same as in the previous figure. In Figure|4]the Hi concentrations in the southern tail 
are labeled SI to S3 and those in the northern tail Nl and N2. Their Hi content is given in Table [3] 
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Fig. 3. Hi spectrum toward the companion IC 4970. 



the 15' beam of the P arkes antenna: 2 0.1 + 0.5 J 
dHorellou & Boothll997l) : 21.7 Jy kms~' dMeyer et 



Jy Km 
al.l2004 



km s 
en 



h-y HIPASS J2016-770 in the HICAT catalog). Since the single- 
dish flux is around 21 Jy km s"', as much as ~ 3.5 Jy km s"' 
could be distributed in an extended area outside the main galax- 
ies and below the detection limit of our interferometric observa- 
tions. This corresponds to as much as ~ 3 x lO^M© of Hi gas, or 
20% the amount of Hi gas contained in the NGC 6872/IC 4970 
system. We can also use the noise level in our maps to estimate 
the Hi flux in the area of the X-ray trail, which iMachacek et al.l 
(12005 ') define as a 435 x 5^9 rectangular region between the 
NGC 6872 and the large elliptical NGC 6876. This area is about 
50 times larger than our synthesized beam. From the noise level 
given in Table 1, we can estimate the noise equivalent flux den- 
sity in that area as NEFD - yfSOa- ^ 6.4 mJy. The radial ve- 
locity difference between the two galaxies is Av ~ 800 km s"', 
which gives a flux NEFD x Av = 5.1 Jy km s"'. This value 
is clearly too high since it exceeds the difference between the 
interferometric and the single-dish measurements. 

The central region of NGC 6872 is devoid of atomic gas. 
Hi was found neither in absorption toward the central continuum 
source, nor in emission. This is not unusual in spiral galaxies, 
where most of the gas in the center is often in molecular form. 
From obser vations of the CO ( J:1-0) line toward the center of 
NGC 6872. lHorellou & BoothI (Il997h inferred a mass of molec- 
ular hydrogen M(H2) of 9.6 x IO^Mq within the central 45". 

The companion galaxy, IC 4970, is clearly detected in 
Hi (Fig. [3]). It contains ~ 1.3 x 10^ Mq of atomic gas. 
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Fig, 4, Grey-scale map of the Hi distribution and Hi spectra measured over several regions in the tails labeled S 1 , S2 and S3 (southern 
concentrations) and NI and N2 (northern concentrations). The first spectrum was measured in the rectangular area around the 
NGC 6872/IC 4970 system. The solid line refers to the ATCA observations, the dotted line is the single-dish Parkes spectrum of 
iHorellou & BoothI ( 119971) and the dotted-dashed line is the Parkes HIPASS spectrum, which has been shifted by -40 mJy for clarity. 
The ten-pointed star s mark the positions of the centers of NGC 6872 and of IC 4970. The location of the star clusters detected by 
iBastian et al.l (l2005h is indicated. Star clusters less massive than IO^Mq are marked by small symbols, and more massive clusters 
by large symbols. Plus signs and squares correspond to star clusters younger than 100 Myr, and crosses and stars of David to older 
clusters. The squares and stars of David indicate a large extinction, with > 1. Note that the regions SI and S2 at the tip of the 
southern tail contain only young star clusters. 
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Fig. 5. The Hi emission of NGC 6872 overlaid on a grey scale optical image from the Digitized Sky Survey and displayed as channel 
maps running from 4496 km s ' (top left panel) to 5116 km s ' (bottom right panel). Every other channel is shown of the channel 
maps showing Hi emission. North is up and east is left. The levels are 3, 5, 7, 9 mJy/beam. The first contour corresponds to the ~ 3cr 
level. The synthesized beam is shown in the bottom left corner of the first panel. 



Five large Hi concentrations are seen in the outer parts of 
NGC 6872, two in the northern tail and three in the southern one. 
Figure |2] shows that the most northern concentration is spatially 
roughly coincident with blue stellar clusters. This is also where 
peaks in the Ho- distribution were seen (Mihos et al. 1993). A 
significant amount of Hi is also found near the break in the north- 
ern arm, close to the nearby companion IC 4970. In the southern 
tail, the three Hi peaks all occur in regions where star clusters are 
seen in the blue-band image. 

There is a clear overall asymmetry in the Hi distribution: we 
estimate that the north-eastern part contai ns about 1.4 times as 
much atomic gas as the south-western part. lBastian et al.l(l2005h 
estimated a star formation rate of ~ 16.5 Moyr"' in the north- 
eastern tail and about half that amount in the western tail, using 
an empirical relationship between the specifi c U -band luminos- 
ity an d the star formation rate measured by iLarsen & Richtleil 
(I2OO O). In the center, they inferred a star formation rate about 
five times lower than in the north-eastern tail. 

Finally, we show integrated Hi spectra measured in different 
regions of the galaxy, as illustrated in Figure|4] The regions in the 
southern tail have been labeled SI to S3, and those in the north- 
ern tail Nl and N2. The amount of gas contained in each region 
is given in Table[3] Gas at the tip of the southern tail in region S 1 
produces a narrow Hi line (width ~ 100 km s ') centered around 
5050 km s"' . Emission from the other regions produces broader 
spectra. In region N2, which lies near the break in the northern 
arm, there are clearly several velocity components. Region Nl 
near the tip of the northern arm also exhibits a narrow Hi line. 

The global Hi spectrum integrated over the whole 
NGC 6872/IC 4970 system is shown in Fig. H (first spec- 
trum, upper left). The line is very broad (~ 950 km s ' at 
the base), in agree r nent with the Parkes sing le-dish Hi spectra 
jHorellou & BoothI 1 19971 iMever et all 120041) . also shown in 



Fig. |4] Since iHorellou & BoothI d 19971) used the radio conven- 
tion Viad/c = Ay/y for velocities, we converted the values to 
show the spectrum using the optical definition v'opt/c = AAjA 
used in the present paper. The low-velocity peak around 
4400 km/s seen in the single-dish Hi spectra is not seen in the 
ATCA map. 

Also shown in Fig. El is the location of the star clusters de- 
tected bv lBastian et al.l (120051) . We have divided the star clusters 
into different groups, depending on their age (more or less than 
100 Myr), mass (more or less than lO^M©), and extinction (more 
or less than A,. = 1), using the estimates of iBastian et al.l (l2005h 
ba sed on the three-dimensional spectral energy fitting algorithm 
of iBik et"aLl (120031) . Interestingly, clusters in the outer regions 
of the tails (SI, S2, Nl, and the region at the very tip of the 
northern tail) are predominantly young and low-mass and have a 
low extinction (small plus signs). Some old, low extinction mas- 
sive clusters (larges crosses) are found north of the northern tail 
at the border of the Hi distribution. Such old massive clusters 
are also found in the inner region (regions S3 and N3). Several 
young, low-mass, low-extinction clusters (small plus signs) are 
found near the companion, slightly to the east, and it is likely that 
their formation has been triggered in the interaction. Regions S3 
and Nl contain large amounts of Hi gas, and a collection of star 
clusters. The young ones (plus signs and boxes) seem to lie at 
the periphery of the Hi concentrations. There is no systematic 
indication that clusters with higher extinction (boxes) lie closer 
to the peak of the Hi concentration. Star clusters are expected 
to form from denser, molecular gas which does not necessarily 
coincide with the more diffuse atomic gas observed here. Also, 
some clusters may lie in the front part of the Hi concent ration and 
not suffer from much extinction. IBastian et al.l (l2005l) found an 
excellent correlation between the location of the youngest clus- 
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Region 


f C J 

J Jhi"v 

T 1 - 1 

Jy km s 


M(HI) 

Mo 


d 

kpc 


NGC 6872/IC 4970 


17.5 


1.54 X 10'" 




NGC 6872 


16.0 


1.41 X 10'" 




IC 4970 


1.46 


1.3 X 10' 




Nl 


1.25 


1.1 X 10' 


17 


N2 


4.40 


3.9 X 10' 


39 


SI 


2.18 


1.9 X 10' 


53 


S2 


1.03 


0.90 X 10' 


34 


S3 


0.86 


0.75 X 10' 


12 



Table 3. Hi fluxes and corresponding Hi masses. The values for 
NGC 6872 were obtained by subtracting those measured toward 
the companion from the values measured for the pair. The third 
column lists the distance on the sky of the peak of the various 
Hi concentrations labeled in Fig.|4]to the center of NGC 6872. 



te rs (< 10 Myr) and t he distribution of Ha emission in the map 
of lMihoset al](ll993h . 

2.2.2. Hi kinematics 

We now examine the kinematical information that can be gained 
from the Hi maps. Figure |5] shows the velocity channel maps as 
isovelocity contours superimposed on an optical image. Gas in 
the north-eastern tail has a lower velocity along the line-of-sight 
(an heliocentric velocity around 4500 km s"') than the south- 
western tail (around 5100 km s"'). The line-of-sight velocities 
vary smoothly across the galaxy and are characteristic of those 
of a rotating disk. 

Figure |6] shows the mean Hi velocity field (top) and the 
Hi velocity dispersion (bottom). Again, the smooth velocity gra- 
dient across the galaxy appears clearly. The velocity dispersion 
map shows a peak on each side of the galaxy's center. One is 
located on the concave side near the break of the northern arm. 
The other one lies south of the main body of the galaxy where 
the southern tail begins. A region of increased velocity disper- 
sion (~ 75 km s"') coincides with the Nl Hi concentration in the 
northern tail. 

3. A^-body simulations 

In order to narrow down the parameter space we first used a re- 
stricted three-body model. 10 000 test particles were distributed 
in a Miyamoto potential representing a galactic disk. The parti- 
cles were advanced according to the gravitational effect of the 
potential of the host galaxy and that of a perturbing companion 
represented by a rigid Plummer potential moving on a Keplerian 
orbit. A good match was found for a prograde, parabolic, low- 
inclination encounter with a companion five times l ess massive 
than th e primary, in agreement with the findings of Mih os et al] 
(Il993h . Then we carried out self-consistent simulations with 
stars and gas, as described below. 

3. 1 . The model 

The simulations were performed using a particle-mesh algo- 
rithm (FFT3D) that makes use of the properties of Fast Fourier 
Transforms and of the convolution theorem. The gravitational 
potential is calculated at the nodes of a three-dimensional 128 x 
128 x64 cartesian grid, providing a resolution of 1 kpc. Effective 
use is made of the entire grid, rather than of only one eighth 




Eight ABaenEion {TSOOO) 

Fig. 6. Top: Hi velocity field in colors superimposed on the VLT 
blue-band image in contours. Bottom: Hi velocity dispersion 
map in colors superimposed on the VLT blue-band image in con- 
tours. 



in standard Fourier algorit hms, t hrough implementation of a 
method described by James (Il977h . 

The main galaxy is made up of three components: a disk, 
a bulge, and a halo. The halo is rigid and represented by a 
Plummer potential, <S>{r) = GM/(r^ + a^y^^. The bulge contains 
particles initially distributed in a Plummer potential. The disk 
contains two types of particles with different velocity disper- 
sions. The collisionless particles representing the stars are ini- 
tially distributed in a Toomre-Kuzmin disk so that the surface 
density has the radial dependance I.(R) - Mcb/[27r(R^ + b^f^]. 
The vertical density profil e corresponds to that of an isothermal 
sheet p(z) oc l/ch^-(z/zQ) (" Spitzejll^il) . The gaseous compo- 
nent is modeled as an ensemble of clouds, all of the same mass, 
that dissipate energy through inelastic collisions with each other 
Gas clouds approaching each other inside a cell collide and ex- 
change energy every tenth time step (every 10^ year). Initially, 
the gas follows a flat radial distribution and accounts for about 
10% of the mass of the stellar disk. The self-gravity of the gas 
is taken into account. Initially, the particles are assigned circular 
velocities in centrifugal equilibrium with the gravitational po- 
tential, a nd velocity dis persions are given by Toomre's stability 
criterion (lToomrelll964 . 

The companion galaxy is modeled as a rigid body repre- 
sented by a Plummer potential. The asumption of a rigid struc- 
ture for the companion is justified insofar as we are interested 
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Fig. 7. Evolution of the stars (left column) and the gas (right column) in the simulation of NGC 6872/IC 4970. The stars in the bulge 
of the primary galaxy are shown in red. The system is viewed according to the observed geometry (position angle of 66°, incUnation 
of 73°). The time step between two snapshots is 40 Myr, starting 20 Myr after perigalacticon. The position of the companion is 
indicated by a star. 
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Orbit Parabolic 

Mass ratio 5:1 

Main galaxy: Mjisk = 120, Odisk = 3 

mass and scale length Mbuige = 30, flbuige = 0.45 

Mbato = 90, flhalo = 10 

orientation ( 1 5 ° , 3 1 5 ° ) 

Companion: mass and scale length M^. = 48, = I 

Table 4. Parameters of the simulations. The orientation of the 
main galaxy is given by the angle between the orbital plane and 
the galactic plane (15°) and the angle between the x-axis in the 
orbital plane and the projection of the vector normal to the main 
galaxy onto the orbital plane. 



in the effect of its passage on the internal structure of the larger 
galaxy, which is probably little affected by the distortion of the 
companion. 

Setting the units of length and time to 1 kpc and 10^ years, 
with the gravitational constant G set to 1, gives a unit of mass of 
2.22 X 10^ Mq. The time step was set to 10*" years, and particles 
were advanced using the leap-frog algorithm. We used 440 000 
particles to model the disk and the bulge of the primary galaxy 
(25% of the particles being initially attributed to the bulge) and 
60 000 gas cloud particles. The parameters of the simulation pre- 
sented here are hsted in TablelH 

3.2. Results 

Figure [T] displays the evolution of the stars and the gas in the 
primary galaxy as it is perturbed by the companion passing on 
a prograde, parabolic orbit. The gas distribution is initially more 
extended than that of the stars. The time between two snapshots 
is 40 Myr, starting 20 Myr after closest approach. Despite its 
small mass (one fifth of the mass of the primary), the companion 
triggers prominent tidal tails very rapidly. The southern tail on 
the side of the galaxy opposite to the companion is thinner that 
the northern tail, where the perturbation is stronger. In the simu- 
lation the companion, modeled as a rigid Plummer potential and 
displayed by a red star, accretes significant amounts of matter 
(both stars and gas). In the observed system, some diffuse stellar 
emission is seen between the companion and the northern arm of 
the primary, but it is not clear that stars or gas have been dragged 
or have fallen onto the companion. Also, in the simulation the 
northern arm breaks at the location of the companion, whereas 
in the observed system both galaxies appear more clearly sepa- 
rated. Aside from those differences, the model is able to repro- 
duce several observed features, including the extent of the tidal 
tails, the thinness of the southern tail, the central bar, the relative 
position of the two galaxies. The simulated system resembles 
the observed one most around 130 Myr after closest approach, 
as shown in Fig. [8] For that epoch, we also display the velocities 
of the gas in the simulations (over the whole galaxy in Fig.|9]and 
in the inner part in Fig. [TOl i. The aiTows indicate the direction 
of the velocity of the gas on the plane of the sky, whereas the 
colors are related to the velocities along the line-of-sight. The 
line-of-sight velocities are in good agreement with the observed 
velocity field. Interestingly, the simulation shows that gas in the 
outer part of the northern tail moves outward, whereas gas in the 
inner part moves inward toward the companion. 

The results of our simulations generally agree with those 
oflMihos etal. (1993). The gas is modeled in a similar way, as 
sticky particles, although in our simulation all gas clouds keep 
the same mass and do not evolve by coalescence or fragmenta- 



tion. Star formation is not included in our model either The ini- 
tial conditions are slightly different: in our simulation, the orbital 
plane of the companion is slightly less inclined on the galac- 
tic plane (15° instead of 25°); their simulated halo is four times 
more massive than the stellar disk, whereas the spherical com- 
ponents of our model are less massive (the added mass of the 
bulge and of the halo is the same as that of the disk). This is why 
the tidal tails are slightly longer and thinner in our simulation. 
We find a blueshif t of the companion relative to the main galaxy 
(-192 km s '), as lMihos et all (1 19931) who found -130 km s '. 
The true velocity difference between the two galaxies is un- 
certain since the published observat ions co v er a l arge range 
of values: For the companion. G reen et alj (Il988b measured 
4759 kms-\ while RC3 give 4727 km s''. For NGC 6872, 
the follow ing values are quoted: 4626 ("Green et al."l988), 4688 
(iMever e t al. 2004, Hi spectrum), 4717 (Horellou & Booth 199^ 
Hi spectrum), 4818 (RC3); 4892 km s ' (IHorellou & BooS 
1997, CO spec trum). For compari s on we converted the veloc- 
ities quoted by iHorellou & Booth] ( Il997h in the radio conven- 
tion Av/ V to the more commonly used optical convention AA/A. 
The broad range of velocities across NGC 6872 (more than 
900 km s"' at the base of the single-dish Hi spectrum) makes 
it difficult to accurately measure the systemic velocity. Taking 
the value from the CO spectrum of NGC 6872, which traces 
the molecular gas in the central 45", and the RC3 value for the 
companion gives a blueshift for the companion of -165 km s"', 
roughly consistent with the simulations. 

Within our model, the various parameters are well con- 
strained by the characteristic morphology of the observed sys- 
tem and we did not find a significantly different set of parameters 
that would give an equally good match to the observations. This, 
however, could change with the addition of more parameters in 
the model (e.g. a halo spin or a different density profile, or an 
additional gravitational perturbation due, for instance, to a third 
galaxy). 

4. Discussion 

Machac ek et al.l ( l2005h discuss several possibilities to explain 
the diffuse X-ray emission seen between NGC 6872 and the 
large elliptical galaxy NGC 6876. 

- As they argue, it is unlikely that the hot gas has been ex- 
tracted from the elliptical galaxy by tidal interactions because 
the mass of the gas in the trail is more than three times larger 
than that contained in the elliptical. Had tidal interactions been 
at work, their effect would have been seen also in the distribution 
of the stars in NGC 6876. The isophotes of NGC 6876 are reg- 
ular; Hubble Space Telescope observations do, however, reveal 
a slight central depression, which is not due to dust absorption. 
This has been interpreted as the possible signature of the end 
product of the merging of t wo gas-free stellar system, each har- 
l5oring a massive black hole (iLauer et al.ll2002l) . We can estimate 
the dynamical mass of N GC 6876 from the observed velocity 
dispersion of 231 km s^' (iDavies et al. l [T987h using the relation 
(M) = (5.0 ± 0.1) X ReCrlliLG), which is the best-fitting virial 
relation found for a sample of 2 5 elliptical and le nticular galax- 
ies from the SAURON sample dCappellari et al.ll2006h . Taking 
the effective radius Re - 25.48 kpc as the half-light radius in 
the blue-band given in NASA/IPAC Extragalactic Database, we 
obtain Mvk = 1.57 x lO'^M©. The projected distance on the sky 
between NGC 6872 and NGC 6876 is ~ 142 kpc, and the radial 
velocity difference is about 800 km s"' . This gives a timescale of 
180 Myr, assuming comparable contributions from the velocity 
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Fig. 8. Distribution the stars (top figure) and the gas (bottom fig- 
ure) in the simulation 130 Myr after perigalacticon, when the 
simulated system resembles the observed one most. The posi- 
tion of the companion is indicated by a star. 
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Fig. 10. Gas in the inner part of the simulated galaxy 130 Myr 
after perigalacticon. The length of the arrows is proportional to 
the velocity on the plane of the sky. The colors indicate the ve- 
locities along the line-of-sight, using the color code as in the 
previous figure. 
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Fig. 9. Gas in the simulation 130 Myr after perigalacticon. The 
arrows indicate the direction of the velocity of the gas clouds 
projected on the plane of the sky. The colors indicate the ve- 
locity range along the line-of-sight, vios: larger than 200 km s"' 
in red, between 100 and 200 km s"' in orange, between and 
100 km s"' in yellow, between -100 and 100 km s"' in green, 
between -200 and -100 km s ' in blue, less than -200 km s ' 
in black. 



in the plane of the sky and the separation along the line-of-sight. 
If NGC 6872 had passed near NGC 6876, a ~ lO'^Mg galaxy 
around 180 Myr ago, it is likely that such an interaction would 
have left a trace in the Hi distribution, which is not observed. 

- The X-ray trail could be due to ram pressure stripping from 
NGC 6872 as the spiral galaxy moves supersonically though 
the hot gas within the Pavo group (Machacek et al. 2005). Ram 
pressure strip ping is known to aff ect the Hi distribution in some 
galaxies (e.g. lVollmer et al.ll2004l) . However, we do not see any 



sign of that effect in the observed Hi distribution of NGC 6872. 
Therefore, we cannot infer any information about the motion of 
NGC 6872 with respect to the surrounding gas in the Pavo group. 

- Other mechanisms involving an interaction between the 
interstellar medium and the intergalactic me dium could b e 
at work, such as turbulent viscous stripping jNulsenl |1982|) . 
where gas at the interface between the galaxy and the inter- 
galactic medium could be heated. At the resolution of our 
Hi observations, we do not see any evidence for that effect ei- 
ther. 

- Accretion into a Bondi-Hoyle trail is the most favored hy- 
pothesis. Machacek et al. (2005) have calculated the e xtent of 
the trail assuming an NFW profile iNavarro et a n (ll996h for the 
dark matter distribution in NGC 6872 with a concentration pa- 
rameter c = 15 and extending until a radius ca, where a is the 
inner radius. An extended halo, with an inner radius a ^ 20 kpc, 
is required to reproduce the observed length of the trail. We have 
performed the same calculation, using a different profile for the 
halo, namely the one corresponding to the Plummer potential 
that we have used in the simulation with a scale length of 10 kpc. 
As illustrated in Fig.[TT] for a Plummer profile the downstream 
profile of density contrast peaks at a larger distance from the 
moving halo, and decreases more softly. At a downstream dis- 
tance z = -lOfl, the density contrast is slightly higher for the 
Plummer sphere than for the NFW profile. Although the density 
contrast along the trail is, in principle, sensitive to the density 
profile of the moving halo, this is difficult to use in practice. 
This calculation is simplified, based on linearized flow equa- 
tions, which applies only for small density perturbations and 
Machacek et al. (2005) point out that the observed temperature 
ratio between the trail and the intra-cluster gas implies an over- 
density larger than one. 
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Fig. 11. Density contrast 5plp downstream of a Plummer po- 
tential (solid lines) and an NFW potential (dotted lines) moving 
supersonically through a medium with different Mach numbers. 

More sensitive X-ray observations would help elucidate the 
nature of the X-ray trail. As discussed above, gravitational in- 
teraction with the large elliptical is unlikely to have produced 
such an X-ray trail without any Hi counterpart. More detailed 
simulations of the evolution of a spiral galaxy in the gravita- 
tional potential of a group, interacting with the hot intra-cluster 
medium and with both a massive elliptical and a smaller galaxy 
are required to shed more light on the relative importance of the 
different mechanisms. This is beyond the scope of this paper. 

5. Summary and conclusion 

In this paper we have presented interferometric Hi observations 
of the interacting system NCG 6872/IC 4970 in the southern 
group of galaxies Pavo obtained with the ATCA, and A^-body 
numerical simulations of the evolution of the system. 

1. We detected about 1.54 x IO^'Mq of Hi gas toward the 
NGC 6872/IC 4970 pair, out of which about 1.3 x lO'^Mg 
is associated with the companion IC 4970. The gas in 
NGC 6872 is distributed in an extended rotating disk. 
Massive concentrations of gas (~ IO^Mq) are found at the 
tip of both tidal tails and towards the break seen in the op- 
tical northern arm near the companion. An increased ve- 
locity dispersion is observed on each side of the galaxy's 
center, near the break in the northern arm and south of the 
main body where the southern tail starts; another region of 
higher velocity dispersion coincides with the most northern 
Hi concentration. We have compar ed the Hi distr i bution with 
that of stellar clusters observed bv lBastian et alj ( 12005 ^. Star 
clusters in the outer Hi concentrations in the tails are pre- 
dominantly young and low-mass (less than 100 Myr and less 
than IO^Mq) and have a low extinction (A,, < 1). 

2. No Hi gas was found associated with the X-ray trail ob- 
served by iMachacek et al.1 (l2005h between NGC 6872 and 
NGC 6876, the dominant elliptical galaxy in the Pavo group 
located ~ 8' to the southeast. However, as much as ~ 3 x 
lO^M© of Hi gas could be hidden in extended regions of dif- 
fuse emission below our detection limit. At the sensitivity 
and the resolution of the observations, there is no sign in 
the overall Hi distribution that NGC 6876 has affected the 
evolution of NGC 6872. There is no evidence either of ram 
pressure stripping. 



3. The simulations of a gravitational interaction with the small 
nearby companion IC 4970 on a low-inclination prograde 
passage are able to reproduce most of the observed fea- 
tures of NGC 6872, including the general morphology of the 
galaxy, the inner bar, the extent of the tidal tails and the thin- 
ness of the southern tail. 

4. Our simulations do not model the gas environment within 
the group or the motion of NGC 6872 through the group. 
Without those interactions we cannot address the hydro- 
dynamic processes, i.e. ram pressure stripping, turbulent 
viscous stripping, or Bondi-Hoyle accretion, discussed in 
Section 4. However, together with the Hi maps, our simula- 
tions provide important constraints on viewing angle, inter- 
nal velocities, mass and gas distributions for input into future 
hydrodynamic simulations. 
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